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摘 要 : 使 用 耦合 化 学 模块 的 高 分 辨 率 中 尺度 数值 模式 WRF-Chem3.4, 结 合 近 地 层 观测 资料 评估 YSU .MYJ、 
QNSE MYNN2.5 和 BouLac 共 5 种 边界 层 参数 化 方案 对 2007 年 3 月 27 日 西北 地 区 一 次 沙 侍 天 气 过 程 模拟 效果 的 影 
响 ,结果 显示 $ 种 边界 层 参 数 化 方案 均 可 模拟 出 此 次 沙 侍 天 气 的 发 展演 变 过 程 ,其 中 YSU 和 BouLac 方 案 模拟 出 相 


对 较 高 的 地 表 摩 擦 速度 、10 m 风速 ,2 m 温 度 和 地 面 PMw 浓 度 以 及 相对 较 低 的 2 mm 相对 湿度 ,从 而 模拟 的 地 表 沙 尘 
天 气 过 程 较 强 ,MYJ、QNSE 和 MYNN2.5 方 案 模拟 的 地 表 沙 侍 天 气 则 相对 较 弱 ,这 表明 不 同 边界 层 参 数 化 方案 通过 
摩擦 速度 的 不 同 模拟 效果 对 沙 侍 排 放 通 量 和 PM 浓度 的 模拟 有 重要 影响 , 较 大 的 摩擦 速度 会 使 起 沙 参数 化 方案 计 
算 的 沙 侍 排 放 通 量 和 PM 浓度 更 高 ,加 之 午后 近 地 层 的 强风 、 高 温和 低 湿 特 征 对 沙 侍 天 气 的 增强 作用 ,使 得 BouLac 


方案 模拟 的 沙 侍 天 气 最 强 ,而 QNSE 方 案 的 模拟 结果 最 弱 ; 利 用 民 勤 站 观测 资料 对 5 种 边界 层 参数 化 方案 模拟 结果 


的 统计 分 析 表 明 ,不 同方 案 对 民 勤 站 沙尘暴 前 后 有 关 气 象 要 素 的 模拟 效果 存在 一 定 的 差异 ,其 中 QNSE 方 案 对 PM 
浓度 的 模拟 效果 最 好 ,BouLac 方 案 对 10 m 风速 的 模拟 效果 最 好 ,YSU 方 案 对 2 m 温 度 和 2 m 相对 湿度 的 模拟 效果 最 


好 ,整体 而 言 ,YSU 方 案 对 民 勤 站 近 地 层 气象 要 素 的 模拟 有 一 定 的 优势 ,QNSE 方 案 的 模拟 结果 相对 最 差 。 


关键 词 : WRF-Chem; WE; 边界 层 参数 化 ; 数值 模拟 


大 气 边 界 层 (Planetary Boundary Layer, PBL) 是 

直接 受 地 表 摩 擦 、 热 力 和 莱 发 过 程 显 著 影响 的 大 

^UZ ,高 度 可 达 数 公里 "。 作 为 自由 大 气 和 地 表 之 

间 的 缓冲 地 带 , 边 界 层 内 由 于 满 流 运动 引起 的 动 
E .热量 和 水 汽 交 换 对 天 气 和 和 气候 有 重要 影响 。 

由 于 利用 数学 方法 描述 油 流 输送 过 程 存在 明 

显 的 困难 ,在 数值 模式 中 通常 以 参数 化 的 形式 来 描 


如 陆 面 过 程 参数 化 、 地 形 和 模式 分 辩 率 等 综合 影 
响 , 局 地 方案 在 模拟 近 地 层 气象 要 素 时 也 有 一 定 的 
优势 ,例如 MYNN2.5 (Mellor-Yamada-Nakanisfii and 
Niino 2.5) 方 案 能 够 较 好 的 模拟 恒 河 流域 一 次 干旱 
过 程 中 。 此 外 ,已 有 研究 表明 边界 层 参数 化 方案 的 
适用 性 因 地 区 差异 、 天 气 过 程 类 型 不 同 存 在 明显 差 
异 。 已 有 学 者 聚焦 局 地 复杂 天 气 过 程 如 台风 、 暴 


述 边 界 层 中 有 关 动 力 及 热力 结构 的 特征 变化 。 不 
同 的 边界 层 参 数 化 方案 因 其 闭合 方法 .闭合 阶 数 等 
方面 的 差异 ,其 数值 模拟 效果 也 存在 一 定 差异 。 相 
关 人 研究 指出 ,与 局 地 方案 相 比 , 非 局 地 方案 通常 具 
有 更 好 的 模拟 效果 ,能够 模拟 出 更 强 的 油 流 混合 
作用 和 更 为 深厚 的 对 流 边界 层 " ,其 中 ACM2 
(Asymmetric Convective Model Version 2) 方 案 是 模拟 
干旱 地 区 大 气 边界 层 结 构 和 近 地 层 气象 要 素 的 较 优 
选择 ” ,可 为 空气 质量 模式 提供 合理 的 气象 场 ” ， 
从 而 更 好 地 模拟 污染 物 浓度 … ,但 考虑 到 其 他 因素 
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雨 .平流 筋 和 龙卷风 等 ,探寻 其 对 不 同 天 气 过 程 的 
影响 。 这 些 研 究 指 出 :ACM2 方 案 和 MYJ(Mellor- 
Yamada-Janjic ) 方 案 对 台风 过 程 的 模拟 有 一 定 的 优 
34/075: YSU (Yonsei University ) 方 案 对 平流 雾 过 程 
中 风速 的 模拟 效果 最 佳 ,QNSE (Quasi-Normal Scale 
Elimination) 方 案 对 平流 雾 生成 的 时 间 地 点 预报 较 
准 “; 非 局 地 方案 低估 了 热力 相关 的 变量 ,局 地 方 
案 可 模拟 出 较 强 的 浮力 ,因而 适合 龙卷风 过 程 的 预 
测 " ;YSU 方 案 适 用 于 四 川 命 地 较 大 量 级 降水 的 模 
拟 "“。 男 有 人 研究 表明 ,同一 参数 化 方案 对 不 同 气 象 
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要 素 的 模拟 效果 也 存在 差异 ””。 以 上 人 研究 的 开展 
对 区 域 数值 天 气 预报 准确 性 的 提高 有 重要 意义 。 
沙 洁 天 气 主要 发 生 在 春季 影响 我 国 西北 地 区 
的 一 种 严重 气象 灾害 ,是 大 风 与 下 垫 面相 互 作用 的 
产物 。 西 北 地 区 沙漠 广 布 ,可 为 此 类 天 气 的 形成 提 
供 充足 的 沙 源 ,春季 大 气 环流 调整 阶段 频繁 南下 的 
冷 空 气 和 西北 地 区 复杂 地 形 的 作用 ,尤其 是 河西 走 
廊 的 狭 管 效应 ,有 利于 大 风 的 形成 和 起 沙 过 程 的 发 
展 ;此 外 ,午后 形成 的 近 地 层 不 稳定 层 结 能 为 沙 尘 
骏 的 形成 提供 有 利 的 不 稳定 机 制 。 有 关 沙 侍 天 气 
的 数值 模拟 研究 常 聚焦 于 数值 模式 中 起 沙 参数 化 
方案 对 起 沙 过 程 的 合理 描述 , 近 地 层 大 气 状态 、 土 
壤 性 质 和 地 表 植 被 等 因素 均 会 影响 起 沙 通 量 的 模 
拟 效 果 ”""。 但 起 沙 过 程 主要 发 生 在 边界 层 ,伴随 着 
风速 . 温 湿 度 、 沙 洁 浓 度 等 要 素 的 显著 变化 ,是 一 种 
有 具有 非 线 性 .多 尺度 、 随 机 性 的 复杂 潮流 输 运 过 程 
,此 外 ,边界 层 满 流 相干 结构 也 与 泥 沙 起 动 与 输 
送 密 切 相 关 ”。 总 而 言 之 ,不 同 边 界 层 参 数 化 方案 
对 满 流 作用 的 模拟 差异 ,会 导致 起 沙 通 量 的 模拟 效 
果 产 生 差异 ,因此 边界 层 参 数 化 方案 的 合理 选取 对 
于 沙 尘 天 气 发 展 过 程 中 沙 侍 浓度 的 模式 预报 有 重 
影响 。 考 虑 到 沙 侍 天 气 的 模拟 与 边界 层 过 程 的 
准确 模拟 密切 相关 ,本文 针对 一 次 西北 地 区 沙 尘 天 
气 过 程 , 主 要 分 析 WRF(Weather Research and Fore- 
casting Model) 模 式 中 有 关 边 界 层 参 数 化 等 对 近 地 
层 气 象 要 素 以 及 PMio 浓 度 模拟 效果 的 影响 。 


1 试验 设计 和 参数 化 方案 简介 


1.1 模式 设置 及 试验 设计 
选择 2007 年 3 月 27 日 西北 地 区 一 次 沙 人 尘 天 气 


90?E 95°E 100?E 105?E TIOE 115?E 


图 1 以 民 勤 站 为 中 心 的 三 重 骨 套 网 格 


99°E 


过 程 为 研究 对 象 ,利用 耦合 了 化 学 模块 的 高 分 辩 率 
中 尺度 数值 模式 WRE- Chem3.4(Weather Research 
and Forecasting Model with Chemistry) 版 本 开展 数值 
模拟 ,重点 关注 非 局 地 闭合 的 YSU 及 局 地 闭合 的 
MYJ、QNSE .MYNN2.5 和 BouLac 共 5 种 边界 层 参数 
化 方案 对 沙 侍 天 气 过 程 中 起 沙 效果 的 影响 。 试 验 
模拟 区 域 中 心 点 设 在 民 勤 站 (103.08°E,38.63°N)， 
采用 三 重 双向 垦 套 网 格 (D01,D02,D03), 模 拟 区 域 
及 D03 地 形 如 图 1 所 示 。 初 始 场 和 边界 条 件 采用 美 
国 气象 环境 预报 中 心 (National Centers for Environ- 
mental Prediction, NCEP) 和 美国 国家 大 气 研究 中 心 
(National Center for Atmospheric Research , NCAR ) 联 
合 制作 的 每 6 一 次 的 1°xX1°FNL(Final) 客 观 分 析 资 
料及 NCEP 全球 地 形 资料 。 上 有 具体 的 模式 设置 及 试验 
设计 在 表 1 中 给 出 。 化 学 模块 的 积分 步 长 为 10 
min, 化 学 方案 选取 仅 考虑 沙 侍 的 GOCART(Global 
Ozone Chemistry Aerosol Radiation and Transport ) 气 
溶胶 方案 ,并 选取 GOCART- AFWA (GOCART with 
Air Force Weather Agency) 方 案 作 为 起 沙 参数 化 方 
案 , 同 时 仅 考 虑 干 沉降 过 程 , 暂 不 考虑 湿 沉 降 、 光 化 
学 过 程 .生物 气 溶 胶 海盐 气 溶胶 . 湿 清 除 及 气 溶 胶 
对 辐射 的 反馈 等 过 程 。 模 式 模 拟 时 D01 模 拟 域 的 时 
间 步 长 为 60 s, 模 拟 时 段 为 2007 年 3 月 27 H 00:00 
UTC 至 28 H 00:00 UTC ,文中 重点 分 析 沙 侍 天 气 过 
程 主要 出 现在 民 勤 站 的 时 段 , 即 27 H 06:00 UTC 至 
15:00 UTC。 
1.2 边界 层 参数 化 方案 

本 次 研究 选取 的 5 种 边界 层 参 数 化 方案 中 ， 
YSU 方 案 “ 采 用 1 阶 闭合 非 局 地 K 理 论 , 采 用 相似 
理论 计算 注 流 扩散 系数 ,适合 模拟 强风 状态 下 的 边 
界 层 结构 ;MYJ、QNSE MYNN2.5 和 BouLac 方 案 均 
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区 域 (a) 及 D03 区 域 的 地 形 高 度 (b) 


Fig. 1 Regional model domains centered at Minqin(a) and spatial distribution of terrain height in D03(b) 
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表 1 模式 设置 及 试验 设计 


Tab.1 Model configuration and experiment settings 


模拟 域 DOI D02 D03 
模拟 时 段 2007 年 3 月 27 日 00:00 UTC-28 H 00:00 UTC 
模拟 中 心 103.08°E ,38.63°N 103.08°E ,38.63°N 103.08°E ,38.63°N 
水 平分 辨 率 15 km 5 km 1.67 km 
格 点 维 数 151x111 271x211 391x331 
垂直 层 数 30 层 30 层 30 层 
顶部 气压 50 hPa 50 hPa 50 hPa 
积分 步 长 60s 20s 6.7 s 
输出 时 间 分 辨 率 Ih 5 min 5 min 
微 物理 过 程 WSMS WSMS WSMS 
短波 辐射 RRTM RRTM RRTM 
长 波 辐射 RRTM RRTM RRTM 
具 云 参数 化 Kain-Fritsch - z 
陆 面 过 程 Noah Noah Noah 
边界 层 及 近 地 层 参数 化 YSU(MM5 M-0),MYJ(M-O) ,QNSE(QNSE) .MYNN2.5(MYNN) 和 BouLac(MM5 M-O) 
起 沙 参 数 化 GOCART-AFWA GOCART-AFWA GOCART-AFWA 


为 1.5 阶 局 地 闭合 方案 ,采用 清流 动能 TKE(Turbu- 
lent Kinetic Energy) 计 算 汕 流 扩散 系数 和 边界 层 顶 
高 度 ,能 预报 局 部 的 垂直 混合 强度 和 潮流 动能 。 
MYJ 方 案 55 采 用 MellorYamada 闭 合 模型 ,在 稳定 条 
件 和 能 不 稳定 条 件 的 边界 层 下 误差 小 ,在 对 流 边界 
层 中 误差 较 大 ,适合 于 精细 边界 层 结构 的 研究 ; 
MYNN2.5 7; && ^* yi] Je JA Mellor-Yamada 闭合 模型 
基础 上 改进 的 预报 次 网 格 动能 方案 ;QNSE DRO 
采用 不 同 的 方案 表征 稳定 和 不 稳定 层 结 下 的 满 流 
Jj Bt. x& HET fe E Ss P ga iiU) 2E BS PSU ; BouLac 
Jr 2€ ^ 3&8 H T BE Uti HU JÉ A B 2 dra Dto ERU BEI 
模拟 和 预报 。5 种 方案 的 主要 特征 如 表 2。 
1.3 起 沙 参 数 化 方案 

本 文选 取 的 起 沙 参 数 化 方案 为 COCART-AF- 


WA (Air Force Weather Agency ) 7 2& ,该 方案 中 粒 径 

按 0.2~2 pm(dust1) 、2~3.6 pum(dust2) ,3.6-6 jum(dust3) , 

6-12 um(dust4) , 12-20 jum(dust5) JE £7 4) P5 , 垂直 沙 
通 量 计算 公式 如下: 

F-aQE (1) 

式 中 : acc 10747797 ,为 考虑 土壤 粘土 比例 的 影响 

的 沙 尘 受 击 率 ; 无 为 风蚀 度 指 数 ;O 为 跃 移 水 平 沙 通 
量 , 采 用 如 下 公式 计算 : 

2 
o-cPafis 30 3l Us S Ws, (2) 


Us. Ux 


AP: C 为 可 变动 的 经 验 比例 常量 ; p, 为 空气 密度 ; 
g 为 重力 加 速度 ; u, 为 摩擦 速度 ; u, 为 临界 摩擦 速 
度 。 从 该 计算 公式 中 可 以 看 出 ,起 沙 通 量 与 边界 层 


表 2 边界 层 参数 化 方案 主要 特征 


Tab.2 Key features of five boundary layer schemes 


报 洪流 动能 强度 


边界 层 方案 闭合 方法 及 阶 数 主要 特征 边界 层 顶 
ysy” 非 局 地 1 阶 增加 反 梯 度 项 及 卷 夹 通 量 项 ,增加 热力 诱导 的 油 流 混合 、 理 查 逊 数 (局 ,不 稳定 层 结 , Ri=0.25， 
` 降低 动力 强迫 的 清流 混合 稳定 ,Ri=0 
Bs Ex 7 反 梯 度 扩散 ,用 Mellor-Yamada 闭 合 模型 表示 表面 层 以 上 滑 流 动能 廓 线 ,TKE 降 至 0.2 ms”? 
MY]? b 1.5 Dj Kiss pene 
fr] aa v Hr eS BE 
不 稳定 层 结 :采用 引入 波动 项 的 MYJ 方 案 清流 动能 闭合 oa aeae 
T EXE y oa Ar MiLa HE SZ , TKE KÆ 0.01 m^ «s? 
QNSE” 局 地 1.5 阶 方案 ; 移 定 层 结 :采用 稳定 条 件 下 改进 的 注 流 谱 闭 合 模式 “ 注 流 可 能 训 线 TREREE OI m?s 
扩散 方案 IRA 
—P — 反 梯 度 扩 散 , 对 气压 相关 项 进行 参数 化 处 理 ,增加 表征 浮 ” 汕 流动 能 廓 线 ,TKE 降 至 10x10-。 
g 力 与 切 变 作用 的 一 组 闭合 参 妆 mès H E E 
Boul ac? 局 地 1.5 阶 能 预报 不 同类 型 下 热 面 潮流 强度 和 具体 位 置 , 并 持续 预 “潮流 动能 廊 线 ,TKE 降 至 0.005 m+s? 
ouLac 可 地 1. 
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中 的 摩擦 速度 u, 及 其 临界 值 密切 相关 ,这 意味 着 不 。 度 ,1 m、2m 和 20m 共 3 层 的 风向 ;采样 间隔 为 300 s 
同 边界 层 参数 化 方案 计算 的 对 起 沙 通 量 的 模拟 。 ”的 PM 浓度 、 环 境 气 压 和 温度 。 
有 重要 影响 。 2007 年 3 月 27 一 28 日 ,中国 西北 地 区 相继 出 现 
由 于 所 选 方案 划分 的 粒 径 范围 与 PM 并 不 完 ”大 范围 沙 企 天 气 , 民 勤 站 监测 到 一 次 沙尘暴 过 程 。 
全 一 致 ,为 了 得 到 与 观测 PM 浓度 对 比 的 模拟 PM。 ”27 H 06:00 UTC, 河 西 走廊 部 分 测 站 出 现 扬 沙 天 气 ， 
ESYA x Yay 5 Ha Ln y N . PAPA y HEP EE 
We RREAUERS IERI UO AEEUEL PMMEUE, np 27 H 0840 UTC, BIO HARDEE, PMR 3S 
算 公 式 如 下 : 9500 ugem”, BERT 10 m 风 速达 到 8 级 (19.4 mes”), 
(10— 6) 能 见 度 不 足 500 m;27 H 12:00 UTC 沙尘暴 东 移 , 河 
Cu, 7 Cue + Casa + Cas + (2-6) Camis (3) AARRE KZAM DU EYERE , B 
E C ond 十 Cino + C us 十 0.67C on 后 沙尘暴 强度 减弱 并 逐渐 消失 ( 图 2) o E 3] X 的 
ee N Ga an a a a a 
P dens » 冷 锋 过 境 是 影响 此 次 沙 尘 天 气 发 生发 展 的 重要 因 
Ea qr—P j: Hm A yhy HS e » v 、 a 
EEE EE E e Ro ETEO RO E AKDERE 
VAM TL . 形 近 地 层 要 素 的 变化 进行 了 细致 分 析 , 故 不 
2 观测 资料 及 沙 尘 天 气 过 程 层 要 素 的 变化 进行 了 细致 分 析 , 故 不 再 
ya MT o 
本 文 将 主要 利用 民 勤 站 近 地 层 观测 资料 来 评 
ae ww 3 模拟 结果 分 析 
佑 不 同 边 界 层 参 数 化 方案 对 一 次 沙 尘 天 气 过 程 的 
影响 。 民 勤 地 处 石 羊 河 最 下 游 ,被 巴 丹 吉林 沙漠 和 31 风 场 变化 分 析 
腾 格 里 沙 江 包 围 , 是 我 国 沙 侍 枯 多 发 区 之 一 。 民 勤 近 地 层 大 风 是 起 沙 过 程 和 沙 尘 天 气 发 展 的 重 
站 (103.08"E,38.63"N) 位 于 河西 走廊 ,海拔 高 度 。 要 因子 。 分析 2007 年 3 月 27 日 5 种 边界 层 参数 化 
1376.5 m, 其 观测 资料 包括 :采样 间隔 为 10s 的 1m、 方案 模拟 的 10 m 风 场 发 现 ,06:00 UTC 模 拟 的 10 m 
2m4m、10m 和 20m 共 5 层 的 风速 温度 和 相对 湿 。 风速 在 民 勤 以 西 出 现 了 风速 大 于 12 ms 的 强风 区 
(a) 06:00 UTC (b) 09:00 UTC 
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* 强 沙 尘 暴 。 民 勤 站 


图 2 2007 年 3 月 27 日 06:00 UTC .09:00 UTC 12:00 UTC 及 15:00 UTC 地 面 监测 的 沙 尘 天 气 
Fig.2 Surface observations from 06:00 UTC to 15:00 UTC on March 27, 2007 with 3 hours interval 
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(图 略 ), 到 09:00 UTC 模 拟 的 10 m 风速 强 风 区 向 东 
扩展 增强 ,并 在 民 勤 以 北 、 以 南 及 以 东 地 区 分 别 形 
成 了 3 个 强风 区 (A、B 和 C), 最 大 风速 大 于 16 ms, 
民 勤 站 位 于 强风 区 B 内 大 值 区 附近 ,5 种 模拟 结 
中 ,YSU 方 案 模 拟 的 风速 中 心 分 布 较为 分 散 ,MYJ 和 
QNSE 方案 模拟 的 风速 中 心 分 布 特征 相似 ， 
MYNN2.5 与 BouLac 方 案 模 拟 的 强风 区 分 布 紧密 ,位 
置 相似 ,但 BouLac 方 案 模拟 的 B 和 C 合 并 ,形成 了 
一 条 东西 走向 的 强风 带 (图 3a)。12:00 UTC,5 种 方 
案 模拟 的 强风 区 向 东南 方向 移动 , 民 勤 已 位 于 强风 
区 西 侧 边 缘 ,其 中 ,YSU 和 BouLac 方 案 模拟 的 强风 
区 持续 发 展 ,MYJ、QNSE 及 MYNN2.5 方 案 模 拟 的 强 
风 区 则 有 所 减弱 ,日 MYNN2.5 方 案 模 拟 的 10 m 风 
速 值 最 低 ( 图 3b)。15:00 UTC,10 m 风速 大 值 区 随 
沙 侍 区 癌 东 南 移动 减弱 ,但 YSU 方 案 和 BouLac 方 案 
模拟 的 10 m 风速 明显 较 大 (图 略 )。 以 上 分 析 表 明 ， 
BouLac 方 案 模拟 的 强风 区 持续 时 间 久 ,YSU 方案 次 
之 ,而 MYJ、QNSE 及 MYNN2.5 方 案 模拟 的 强风 区 持 
续 时 间 短 ,并 在 沙 侍 天 气 发 展 旺 盛 阶段 减弱 ,其 中 
MYNN2.5 方 案 的 减弱 程度 最 为 明显 。 

此 外 ,分 析 民 勤 站 上 空 水 平 风 场 垂直 分 布 的 时 
间 演 变 (图 4) 发 现 ,5 种 方案 均 模 拟 出 了 沙尘暴 发 展 
阶段 高 空 动 量 下 传 造成 的 低层 风速 增加 ,在 09:00 
UTC ,5 种 方案 模拟 的 近 地 层 风速 均 达到 了 20 m-s”, 


(al) YSU 09:00 UTC (a2) MYJ 09:00 UTC 
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Ez: ud 3 
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(bl) YSU 12:00 UTC (b2) MYJ 12:00 UTC 


104*E 106*E 104*E 106*E 104^E 


但 YSU MYJ 和 QNSE 方 案 模 拟 的 高 空 风速 较 强 , 风 
速 高 值 区 (>30 m- s) p] F SE fif £ 2] 3000 m, 
MYNN2.5 和 BouLac 方 案 模拟 的 高 空 风速 较 弱 ,风速 
高 值 区 (>30 ms ) 仅 向 下 延伸 至 4400 nm 左右 。 

3.2 摩擦 速度 场 变 化 分 析 

GOCART-AFW A 起 沙 参数 化 方案 中 ,起 沙 通 量 
的 计算 依赖 于 边界 层 中 摩擦 速度 ,的 大 小 ,因此 不 
同 边界 层 参 数 化 方案 计算 的 对 起 沙 通 量 及 PM 
浓度 的 模拟 有 重要 影响 。 分 析 27 日 5 种 边界 层 参 
数 化 方案 模拟 的 平均 摩擦 速度 发 现 , 民 勤 站 在 
06:00 UTC 位 于 摩擦 速度 大 于 0.9 ms ' 的 高 值 区 
内 ,随后 高 值 区 逐渐 东 移 扩大 ,到 09:00 UTC, 民 勤 
及 其 附近 地 区 的 摩擦 速度 升 高 至 1.2 ms 以 上 , 随 
后 高 值 区 逐渐 减弱 并 移 至 民 勤 以 东 地 区 (图 5) ,在 
此 期 间 ,摩擦 速度 高 值 区 与 沙 尘 区 的 演变 有 较 好 的 
对 应 。 

为 细致 分 析 沙 侍 天 气 发 生发 展期 间 ( 初 期 ; 
06:00 一 09:00 UTC, 中 期 :09:00 一 12:00 UTC, 后 期 : 
12:00—15:00 UTC) 不 同 参 数 化 方案 对 摩擦 速度 的 
模拟 差异 ,将 不 同 参 数 化 方案 模拟 的 摩擦 速度 场 与 
5 种 方案 模拟 结果 的 平均 场 作 差 得 到 不 同时 刻 的 距 
平 场 , 再 做 时 间 平 均 后 得 到 不 同时 段 的 平均 距 平 
场 。 分 析 发 现 ,YSU 方 案 模拟 的 摩擦 速度 在 沙 尘 天 
气 初期 低 于 平均 值 , 但 在 沙 侍 天 气 中 后 期 高 于 平均 


(a3) QNSE 09:00 UTC (a4) MYNN2.5 09:00 UTC (a5) BouLac 09:00 UTC 
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图 3 2007 年 3 月 27 日 不 同时 刻 D03 区 域内 地 面 监测 的 沙 侍 天 气 及 $ 种 边界 层 参 数 化 方案 模拟 的 10 m 风 场 


Fig.3 Surface observations and wind fields at 10 m level simulated by different PBL schemes in D03 at 
different times on March 27, 2007 
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4 20074E3 H 27 H 03:00 UTC £ 15:00 UTC D03 区 域内 5 种 边界 层 参 数 化 方案 模拟 的 民 勤 站 上 空 水 平 风 场 
Fig.4 Vertical distribution of horizontal wind field simulated by five different PBL schemes at Mingin station in 
D03 from 03:00 UTC to 15:00 UTC on March 27, 2007 
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(b) 09:00 UTC 


(c) 12:00 UTC 
102°E 104°E 106°E 
1 L 


41°N 


40°N 


39°N 


38°N 


37°N 4 


102°E 104°E 106°E 


EIKI BE/(m-s) 


ll 1:2 L3 14 15 


图 5 200743 H 27 日 不 同时 刻 D03 区 域内 沙 侍 天 气 地 面 监测 及 5 种 边界 层 参数 化 方案 模拟 的 平均 摩擦 速度 
Fig. 5 Surface observations and average friction velocity simulated by five different PBL schemes in 
D03 at different times on March 27, 2007 


fE ; MY] FI QNSE 77 RIRIA B] EE PK E HE TE EI Vb to 
期 均 偏 低 ;MYNN2.5 方 案 模 拟 的 摩擦 速度 在 沙 尘 天 
气 初期 略 高 于 平均 值 ,但 在 中 后 期 逐渐 低 于 平均 
值 ;BouLac 方 案 模 拟 的 摩擦 速度 在 沙 尘 天 气 初期 、 
中 期 和 后 期 均 高 于 平均 值 ; 在 沙 侍 天 气 发 展期 间 平 
均 而 言 ,YSU 方 案 和 BouLac 方 案 模拟 的 摩擦 速度 高 


3.3 2m 温 湿 场 变 化 分 析 

温 湿 场 的 变化 对 沙 侍 天 气 的 发 展 也 有 重要 作 
用 ,午后 地 表 的 高 温 低 湿 特征 能 够 促进 起 沙 过 程 的 
发 展 ,增强 沙 侍 天 气 的 强度 。 利 用 相同 方法 分 析 2 
年 温 湿 平均 距 平 场 发 现 , 沙 侍 天 气 发 生 期 间 ,BouLac 
和 YSU 方 案 模拟 的 2 m 温 度 高 于 均值 .2 m 相对 湿度 


于 均值 ,其 中 BouLac 方 案 模 拟 的 摩擦 速度 明显 偏 
高 ,MYJ 方 案 、MYNN2.5 方 案 和 QNSE 方 案 模 拟 的 摩 
察 速度 低 于 平均 值 (图 6) ,以 上 特征 均 在 沙 侍 天 和 气 
后 期 最 为 明显 。 


低 于 均值 ,其 中 BouLac 方 案 模拟 的 温度 最 高 .相对 
湿度 最 低 , 且 YSU 方 案 模拟 结果 的 高 温 、 低 湿 特 征 
与 黄 文 谍 等 9 对 该 方案 的 评估 结果 一 致 ;MYJ 方 案 、 
MYNN2.5 方 案 和 QNSE 方 案 模拟 的 2 m 温 度 低 于 均 
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IÈ: (a) X 3 h F94 (06:00-09:00 UTC) , (b) 9 3 h F194 (09:00-12:00 UTC), ()29 3 h FI (12:00-15:00 UTC), 
(d)29 9 h F (06:00-15:00 UTC), 
图 6 2007 年 3 月 27 日 D03 区 域内 5 种 边界 层 参数 化 方案 模拟 的 摩擦 速度 平均 距 平 场 


Fig.6 Average differences between simulated friction velocity and average friction velocity in D03 on March 27, 2007 
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值 .2mm 相 对 湿度 高 于 均值 ,其 中 QNSE 方 案 模拟 的 
温度 最 低 、 相 对 湿度 最 高 , 且 5 种 方案 的 模拟 特征 均 
在 沙 尘 天 气 后 期 最 为 明显 (图 略 ) ,这 与 摩擦 速度 模 
拟 结果 的 特征 一 致 。 
3.4 PM 浓度 场 变化 分 析 

图 7 显示 了 2007 年 3 月 27 H 0853 UTC 气象 卫 
星 监 测 的 沙 侍 区 分 布 及 D02 区 域内 5 种 方案 模拟 的 
27 日 08:55 UTC PM 浓度 空间 分 布 。 卫 星 监测 结果 
显示 , 民 勤 以 北 (A) 及 河西 走廊 (B) 地 区 为 沙 尘 区 ; 
模拟 结果 显示 ,5 种 方案 模拟 的 PM 浓度 高 值 区 位 于 
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WEK A 南 侧 和 区 域 C(102°E~106?°E,38°N~39?°N) 
内 , 沙 尘 区 B 大 部 分 区 域内 模拟 的 PM 浓度 大 于 
250 ug'm” ,表明 5 种 边界 层 参数 化 方案 模拟 的 沙 侍 
区 分 布 与 卫星 监测 结果 相 比 有 一 定 的 合理 性 ,但 各 
方案 模拟 的 沙 尘 区 和 A 均 偏 南 , 沙 侍 区 B 内 的 PM 浓 
度 相 对 较 低 ,此 外 模拟 结果 在 云 区 C 还 出 现 一 个 
PM 浓度 的 高 值 区 。 

对 比分 析 不 同时 刻 地 面 观测 的 天 气 现象 和 5 种 
边界 层 参 数 化 方案 模拟 的 地 面 PMw 浓 度 可 以 发 现 ， 
5 种 方案 均 模 拟 出 了 此 次 沙 侍 天 气 初 生 ` 发 展 及 消 


(b) YSU 09:00 UTC 
98°E 100°E 102°E 104°E 106°E 


98°E 100°E 102°E 104°E 106°E 
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(f) BouLac 09:00 UTC 
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图 7 20074E3 H 27 H 08:53 UTC 气象 卫星 监测 的 沙 侍 区 (9 及 D02 区 域内 08:55 UTC 5 种 不 同 边界 
层 方案 模拟 的 地 面 PM 浓 度 (b~f 
Fig. 7 Dust and sand distribution observed by meteorological satellite at 08:53 UTC (a) and PM» concentration 
(b-f) at surface simulated by five different PBL schemes at 08:55 UTC on March 27 2007 in D02 
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亡 过 程 中 PMw 浓 度 的 演变 。 地 面 监测 显示 ,27 日 
06:00 UTC, 沙 侍 天 气 开始 发 展 ,河西 走廊 及 民 勤 以 
西部 分 测 站 出 现 扬 沙 天 气 ( 图 略 ) ,到 09:00 UTC, 沙 
侍 区 进一步 向 东南 扩展 , 民 勤 以 北 、 以 东 共 2 个 测 站 
出 现 强 沙 侍 暴 ,此 时 民 勤 沙尘暴 已 发 展 到 最 强 时 
Ajo 与 监测 结果 相似 ,5 种 方案 模拟 的 PMio 浓 度 高 


3.5 模拟 和 观测 的 民 勤 观测 站 近 地 层 气象 要 素 的 
对 比分 析 

根据 第 2 节 的 分 析 表 明 , 民 勤 站 在 2007 年 3 月 
27 日 监测 到 一 次 沙尘暴 过 程 ,图 9 显示 了 27 日 民 勤 
站 一 些 气象 要 素 的 地 面 观测 值 及 5 种 边界 层 方案 的 
模拟 值 , 各 要 素 观 测 和 模拟 的 均值 及 峰值 如 表 3 所 


值 区 均 向 东南 方向 扩展 增强 ,并 在 10 m 风 速 强风 区 
A、B 和 C 分 别 形成 了 3 块 沙 尘 区 ,其 强度 和 分 布 特 
征 与 强风 区 基本 一 致 (图 8a)。 此 外 ,5 种 方案 均 模 
拟 出 了 民 勤 站 附近 的 PM 浓度 高 值 中 心 ,其 中 
YSU MYNN2.5 fill BouLac 方案 模拟 的 PM 浓度 较 
高 , 且 BouLac 方案 模 拟 的 PM 浓度 最 大 值 超过 了 
13000 hg'm “(图 8bp)。 进 一 步 分 析 民 勤 站 及 其 附近 
区 域内 小 侍 粒 子 的 三 维 分 布 (图 8c) 发 现 ,BouLac 方 
案 模 拟 的 沙 侍 中 心 (103.1*E,38.6"5N ) 上 空 有 较 强 的 
上 升 气流 , 沙 侍 柱 向 上 延伸 至 4000 m 高 度 以 上 ， 
PM 浓度 高 值 区 (>8000 pug m?) [6] _E fiEffi Æ 2000 m 
高 度 左右 ;YSU 和 MYNN2.5 方 案 模 拟 的 PMw 浓 度 高 
值 区 主要 位 于 低层 ,垂直 伸展 高 度 相对 较 低 ; MYJ 
和 QNSE 方 案 模拟 的 PMio 浓 度 垂直 分 布 特征 相似 ， 
沙 尘 垂直 伸展 高 度 最 低 。12:00 UTC, 沙 尘 天 气 达 到 
最 强 ,河套 以 西部 分 测 站 出 现 沙 尘 暴 及 强 沙 尘 暴 天 
气 ;5 种 方案 模拟 的 PM 浓度 高 值 区 东 移 至 河套 西 
部 并 合并 增强 ,其 中 ,YSU 和 BouLac 方 案 模 拟 的 强 
沙 区 (>10000 pgm ) 分 布 范 围 明 显 较 广 、 强 度 较 强 
(图 8d)。 此 后 ,地 面 监测 的 沙 尘 区 东 移 减弱 ,5 种 方 
案 模 拟 的 PM 浓 度 在 河套 地 区 东部 降 至 最 低 ,其 中 
MYJ 和 QNSE 方 案 模拟 的 PMio 浓 度 降幅 最 大 ,Bou- 
Lac 方 案 降 幅 最 小 (图 略 )。 以 上 分 析 表 明 ,YSU 和 
BouLac 方 案 模拟 的 沙 尘 区 呈现 分 布 范围 广 .持续 时 
间 久 、PMw 浓 度 高 的 特征 ,特别 是 在 沙 尘 天 气 发 展 旺 
威 时 刻 ( 约 12:00 UTC) ,与 地 面 监测 结果 有 较 好 的 
一 致 性 ,其 中 BouLac 方 案 模拟 的 沙 尘 区 最 强 ; 
MYNN2.5 方 案 虽 然 在 09:00 UTC 模拟 出 了 民 勤 的 
PMK RE REX ,但 与 其 他 方案 相 比 ,其 在 沙 尘 天 气 
发 展 旺盛 阶段 模拟 的 PM 浓度 较 低 ;MYJ 方 案 和 
QNSE 方 案 模拟 的 沙 侍 区 分 布 范围 小 .持续 时 间 短 ， 
沙 尘 区 内 PM 浓度 在 沙 人 尘 天 气 发 展 旺盛 阶段 较 
低 。 此 外 ,对 比 PMw 浓 度 与 10 m 风速 的 模拟 结果 可 
知 ,二 者 的 变化 特征 有 和 较 好 的 一 致 性 ,PM 浓度 高 值 
区 总 与 强风 区 相伴 , 且 强 风 区 内 风速 越 高 ,PM 浓度 
越 高 ,这 表明 强风 区 是 影响 强 沙 区 形成 的 主要 
因子 。 


示 ,不 同 要 素 的 平均 时 段 不 同 。PM 浓 度 和 10m 风 
速 的 模拟 值 在 沙尘暴 发 生 之 前 与 观测 值 相差 较 大 ， 
重点 关注 沙尘暴 发 生 期 间 二 者 模拟 值 与 观测 值 较 
为 吻合 的 时 段 (07:00 一 12:00 UTC) ,各 方案 模拟 的 
2 m 温 度 和 2 m 相 对 湿度 在 沙尘暴 发 生前 后 较为 接 
近 且 与 观测 值 相差 较 小 ,因此 重点 关注 二 者 的 日 变 
化 , 沙 尘 排放 通 量 和 地 表 摩 擦 速度 的 平均 时 段 为 沙 
尘 排 放 通 量 的 产生 时 段 (03:00 一 12:00 UTC ) WMI 
的 10 mm 风速 从 05:30 UTC 左右 开始 迅速 上 升 ,并 在 
09:00 UTC 左右 达到 峰值 (15.0 m* s) ,5 种 方案 模拟 
HB 10 m 风速 也 均 在 09:00 UTC 左右 达到 峰值 ,但 在 
沙尘暴 前 期 模拟 值 明显 较 高 (图 9a)。2 m 气 温和 2 
m 相对 湿度 的 观测 结果 与 5 种 方案 的 模拟 结果 基本 
一 致 ,但 14:00 UTC 之 后 ,二 者 的 模拟 值 与 观测 值 出 
现 明显 偏差 ,2 mm 气温 的 模拟 值 较 观 测 值 偏 低 约 
4% ,这 可 能 是 模拟 时 没有 考虑 沙 尘 的 辐射 效应 所 
致 ,2mm 相 对 湿度 的 模拟 值 逐 渐 高 于 观测 值 ,其 中 
YSU 方 案 模拟 的 相对 湿度 与 观测 结果 最 为 接近 ， 
QNSE 方案 模拟 的 相对 湿度 则 较 观 测 结果 高 20% 
(图 9b, 9c)。 在 强风 、 高 温和 低 湿 特征 作用 下 ,观测 
Bj PM, YK BE M 27 H 07:00 UTC 左右 开始 激增 ,在 
08:40 UTC 及 09:40 UTC 先后 达到 9647 ug* m^ 及 
8337 ug m^, BEJATE 12:00 UTC JEJE 4824] 200 jug m^; 
5 种 方案 模拟 的 PM。 浓 度 均 在 09:00 UTC 左右 达到 
峰值 ,上 且 较 观测 结果 提前 约 3 h 开 始 上 升 ,滞后 约 1h 
降 至 500 ug m^ 左右 ;此 外 ,不 同方 案 在 沙尘暴 发 生 
期 间 模拟 的 PM 浓度 大 小 也 存在 差异 ,BouLac 和 
MYNN2.5 方 案 的 模拟 结果 高 于 观测 结果 ,YSU 方案 
的 模拟 结果 与 观测 结果 较为 接近 ,MYJ 和 QNSE 方 
案 的 模拟 结果 低 于 观测 结果 (图 9d ) 

模式 模拟 的 PMw 浓 度 大 小 依赖 于 起 沙 通 量 的 
模拟 ,起 沙 通 量 的 计算 又 与 地 表 摩 擦 速度 密切 相 
关 , 分 析 发 现 ,5 种 方案 模拟 的 沙 尘 排放 通 量 和 地 表 
摩擦 速度 均 在 09:00 UTC 左右 达到 峰值 (图 9e, 9f) , 
1H. 5j MYJ ftl QNSE 7r Z& 4H EG , MYNN2.5 fll BouLac 77 
案 模 拟 的 沙 侍 排放 通 量 和 地 表 摩 擦 速度 的 均值 和 
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TE : (a) J3 09:00 UTC Hh rti PM oT RE FNE rt AUS , (b) 1 BACCHI (a) P 27 EE XR, (e) A dr (b) rr ERKA, Wx20 和 PM 浓 度 垂 直 剖 面 ， 
(gd) 为 12:00 UTC Hh rfi PM o ZEE RU TG XU 
图 8 20074E3 H 27 H 09:00 UTC(a, b, c)ffll 12:00 UTC(d)D03 区 域内 地 面 监测 的 沙 尘 天 气 及 5 种 边界 层 参 数 化 
方案 模拟 的 PM 浓度 及 10 m 风 场 
Fig. 8 Surface observations and PM» concentrations and wind fields at 10 m level simulated by five different PBL 
schemes in D03 at 09:00 UTC(a, b, c) and 12:00 UTC(d) on March 27, 2007 
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图 9 2007 年 3 月 27 日 D03 区 域内 民 勤 站 10 m 风 速 .2m 温 度 .2m 相 对 湿度 和 地 面 PM 浓 度 的 地 面 观 测 和 不 同 边界 层 


方案 模拟 结果 以 及 不 同 边界 层 方 案 模拟 的 沙 尘 排放 通 量 和 摩擦 速度 


Fig.9 Temporal variety of simulated and observed wind speed at 10 m, temperature at 2 m, relative humidity at 2 m and PM, 


concentration at surface, and simulated dust emission, friction velocity at Minqin station in D03 on March 28, 2007 


峰值 较 高 ,YSU 方案 次 之 ( 表 3)。 综 合 分 析 各 方案 YSU、BouLac 和 MYNN2.5 方 案 模拟 的 摩擦 速度 . 沙 
对 地 表 摩 擦 速度 、 沙 侍 排 放 通 量 和 PM 浓度 的 模拟 侍 排 放 通 量 和 PM 浓度 较 高 , MYJ 和 QNSE 方 案 的 
结果 发 现 ,3 种 要 素 呈 现 出 一 致 的 强度 变化 特征 , 即 相关 模拟 结果 最 低 。 以 上 特征 表明 ,在 起 沙 参 数 化 
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表 3 2007 年 3 月 27 日 不 同 气象 要 素 观测 及 模拟 的 均值 和 峰值 


Tab.3 Different observation and simulation of the mean and peak on March 27, 2007 
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气象 要 素 平均 时 段 OBS YSU MYJ QNSE MYNN2.5 BouLac 
均值 。 07:00—12:00 UTC 10.6 10.4 11.2 10.3 11.2 11.3 
10m 风 速 /ms) [t a 15 13.2 15 14 14.3 14.1 
均值 00:00—00:00 UTC 12.6 11.1 10.7 10.3 10.6 11 
2 mili E/'C [t = 21.8 20.5 214 20.7 20.7 212 
均值 。 00:00—00:00 UTC 26 25 29 32 29 27 
2 m 相 对 湿度 /% [t z 43 50 58 65 56 57 
. 均值 07:00—12:00 UTC 4503 4787 3402 3301 5944 6713 
PM 浓度 /(g m?) 
条 值 - 9647 9658 6946 5462 12255 16129 
. mn . 4f 。 03:00 一 12:00 UTC S 694 540 448 1016 861 
VbEAEERIGE SI (ag m^ s) 
条 值 - - 2996 2548 1378 4497 4547 
" ^fi ^ 03:00—12:00 UTC - 0.90 0.86 0.83 0.90 0.90 
摩擦 速度 /(m*s') 
条 值 - = 1.17 1.14 1.08 1.25 1.25 


方案 相同 时 ,不 同 边 界 层 方案 对 摩擦 速度 的 模拟 差 
异 影响 了 沙 尘 排 放 通 量 和 PM 浓度 的 模拟 结果 , 较 
大 的 摩擦 速度 会 使 起 沙 参数 化 方案 计算 出 更 高 的 
沙 侍 排放 通 量 和 较 高 的 PM 浓度 。 

为 定量 评估 5 种 边界 层 参 数 化 方案 对 民 勤 站 近 
地 层 气象 要 素 的 模拟 效果 ,选取 相关 系数 (R) 和 均 
方 根 误差 (RMSE) 对 模拟 结果 进行 检验 ,二 者 的 计 
算 方 法 如 下 : 


(5) 


AP: P, 为 第 i 时 刻 的 模拟 值 ; 0; 为 第 时刻 的 观 
测 值 ; P 为 模拟 的 平均 值 ; 0 为 观测 的 平均 值 ; N 为 
时 间 序 列 的 长 度 。 

5 种 边界 层 参数 化 方案 对 民 勤 站 地 面 PMv 浓 


注 :00:00 一 00:00 UTC 表示 27 H 00:00 UTC 至 28 H 00:00 UTC ; OBS 表示 观测 值 。 


度 10m 风速. 2 m 气 温和 2m 相 对 湿度 模拟 效果 的 
统计 分 析 结 果 可 以 看 出 ( 表 4) ,BouLac 方 案 模拟 的 
10 m 风速 与 观测 值 的 相关 系数 最 高 . 均 方 根 误差 最 
低 ,MYJ 和 QNSE 方 案 模拟 的 10 m 风 速 与 观测 值 的 
相关 系数 最 低 均 方 根 误差 最 高 。 同 时 ,5 种 边界 层 
参数 化 方案 均 能 较 好 的 模拟 出 局 地 天 气 的 日 变化 
特征 ,模拟 的 2 m 温 度 和 2 mm 相对 湿度 与 观测 结果 的 
相关 系数 均 高 达 0.98。 为 统计 分 析 不 同方 案 模 拟 的 
济 流 特征 差异 ,利用 2 m 温 度 及 2 mm 相对 湿度 的 扰动 
值 来 探讨 不 同方 案 的 模拟 效果 ,扰动 序列 是 由 原始 
序列 减 去 利用 观测 序列 拟 合 出 的 日 变化 趋势 得 
出 。 统 计 结 果 中 ,YSU 方 案 模 拟 的 扰动 2m 温 度 和 
扰动 2 m 相对 湿度 与 观测 值 的 相关 系数 最 高 . 均 方 
根 误差 最 低 ,QNSE 方 案 的 相关 系数 最 低 , 均 方 根 误 
差 最 高 ,此 外 ,与 模拟 和 观测 的 原始 序列 相关 性 相 
比 ,模拟 和 观测 的 扰动 序列 相关 性 大 幅 下 降 , 这 显 
然 与 边界 层 满 流 过 程 的 复杂 性 有 关 。5 种 方案 模拟 
的 PM 浓度 与 观测 值 的 相关 系数 则 差异 较 小 ,但 


表 4 不 同 边 再 层 参 数 化 方案 对 民 勤 站 近 地 层 气象 要 素 及 其 扰动 模拟 结果 的 统计 分 析 


Tab.4 Statistical analysis of near-surface meteorological elements for different PBL schemes at Minqin 


10 m 风速 /ms 


2m 温 度 (扰动 2mm 温 度 )/%C 


2 m 相对 湿度 (扰动 2 m 相 对 湿度 )/% PM IK Eg: m?) 


边界 层 方案 

R RMSE R RMSE R RMSE R RMSE 
YSU 0.86 2.33 0.98(0.24) 2.11 0.98(-0.13) 3.26 0.78 1386 
MYJ 0.84 2.85 0.98(0.16) 2.15 0.97(-0.04) 6.57 0.80 1253 
QNSE 0.79 2.79 0.98(0.13) 3.01 0.97(-0.01) 10.44 0.81 1220 
MYNN2.5 0.89 2.21 0.98(0.15) 2.69 0.98(-0.06) 6.31 0.80 1725 
BouLac 0.92 1.93 0.98(0.23) 2.28 0.96(-0.06) 5.11 0.78 2157 


注 :表示 相关 系数 ,RMSE 表示 均 方 根 误 差 ;* 表 示 未 通过 a=0.05 的 显著 性 检验 。 
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YSU 、MYJ fll QNSE 方案 模拟 的 PMw 浓 度 与 观测 值 
的 均 方 根 误差 较 低 ,MYNN2.5 和 BouLac 方 案 模拟 的 
PMi 浓 度 与 观测 值 的 均 方 根 误差 较 高 。 

以 上 分 析 表 明 ,就 民 勤 站 沙尘暴 前 后 地 面 气象 
要 素 的 模拟 而 言 ,不 同方 案 对 同一 气象 要 素 的 模拟 
效果 不 同 , 从 统计 结果 来 看 ,BouLac 方案 对 10m 风 
速 的 模拟 效果 最 好 ,YSU 方案 对 2 m 温 度 和 2 mm 相对 
湿度 的 模拟 效果 最 好 ,QNSE 方案 对 PM 浓度 的 模 
拟 效果 最 好 。 但 综合 分 析 各 方案 对 民 勤 站 4 种 地 面 
气象 要 素 的 模拟 效果 发 现 ,YSU 方 案 对 民 勤 站 地 面 
气象 要 素 的 模拟 有 一 定 的 优势 ,QNSE 方 案 对 这 些 
要 素 的 模拟 效果 则 相对 较 差 。YSU 方 案 对 2m 温 度 
及 相对 湿度 的 模拟 均 呈 明显 优势 , 且 对 10m 风速 的 
模拟 效果 较 好 ,最 终 模 拟 的 PM' 浓 度 值 与 观测 值 的 
误差 小 , 且 在 沙尘暴 期 间 模 拟 的 平均 浓度 和 峰值 浓 
度 与 观测 最 为 接近 ( 表 3)。QNSE 方 案 模拟 的 10 m 
风速 .2 m 温 度 和 2 m 相对 湿度 与 观测 值 比 相关 性 
低 、 误 差 大 , 且 其 模拟 的 平均 和 峰值 PM, 浓 度 与 观 
测 结果 相差 较 大 。 

本 文 针 对 西北 地 区 一 次 沙 侍 天 气 过 程 初步 评 
估 了 不 同方 案 模拟 沙 尘 天 气 的 主要 特征 和 差异 ,但 
并 未 对 差异 产生 的 具体 原因 做 过 多 研究 。 由 于 地 
表土 壤 的 风蚀 起 沙 过 程 是 边界 层 大 气 与 地 面 沙 尘 
粒子 相互 作用 的 一 种 动力 过 程 ,与 边界 层 内 复杂 的 
满 流 输 运 过 程 密 不 可 分 ,因此 进一步 深入 研究 不 同 
边界 层 参数 化 方案 模拟 结果 差异 的 产生 原因 和 边 
界 层 过 程 对 沙 侍 输送 机 制 的 影响 可 作为 下 一 步 研 
究 方 向 ,这 对 于 提高 沙 尘 天 气 的 模拟 和 预报 有 重要 


| 
AA o 


4 结论 


本 文 使 用 耦合 化 学 模块 的 高 分 辨 率 中 尺度 数 
值 模式 又 RF-Chem3.4, 结 合 观测 资料 评估 不 同 边界 
层 参数 化 方案 对 一 次 沙 尘 天 气 过 程 的 模拟 效果 ,并 
对 比分 析 了 5 种 不 同 边界 层 参数 化 方案 模拟 的 近 地 
层 气象 要 素 特 征 的 有 关 异 同 ,得 出 以 下 结论 : 

(1) 选取 5 种 不 同 边 界 层 方案 的 WRF-Chem3.4 
模式 均 在 中 国 西北 地 区 模拟 出 了 此 次 沙 侍 天 气 的 
初生 、 爆 发 及 减弱 过 程 ,与 MYJ 方 案 、QNSE 方 案 及 
MNN2.5 FRH Eb , YSU 方案 和 BouLac 方 案 的 模拟 
结果 沙 侍 区 内 表现 为 强风 速 和 高 PM 浓度 的 特征 ， 
且 能 模拟 出 沙 人 尘 天 气 发 展 旺盛 阶 段 更 高 的 2 m 温 度 


及 更 低 的 2 m 湿 度 等 特征 ,因此 这 两 种 方案 模拟 出 
的 沙 侍 天 气相 对 较 强 , 这 与 地 面 监测 结果 有 较 好 的 
— SB lE; QNSE 方案 模拟 出 较 其 他 几 种 边界 层 参数 
化 方案 更 低 的 10 m 风速 .PM 浓度 和 2 m 温 度 以 及 
更 高 的 2 mm 湿度 。 

(2) 地 面 风速 和 PM 浓度 的 模拟 结果 表明 , 强 
风 区 是 导致 PM 浓度 高 值 区 形成 的 主要 因子 ,强风 
区 内 风速 越 高 PM 浓度 越 高 ;地 表 摩 擦 速度 、 沙 侍 
排放 通 量 和 PM 浓度 的 模拟 结果 表明 ,摩擦 速度 是 
影响 PM 浓度 高 值 区 形成 的 重要 因子 ,因此 不 同 边 
界 层 方案 可 通过 对 摩擦 速度 的 不 同 模拟 效果 影响 
沙 侍 排放 通 量 和 PM 浓度 , 较 大 的 摩擦 速度 会 使 起 
沙 参 数 化 方案 计算 更 高 的 沙 侍 排放 通 量 ,这 使 YSU 
方案 和 BouLac 方 案 模 拟 出 较 高 的 PM 浓度 。 

(3) 结合 民 勤 站 观测 资料 进行 统计 分 析 得 出 ， 
不 同方 案 对 民 勤 站 沙尘暴 前 后 同一 地 面 气象 要 素 
的 模拟 效果 不 同 ,YSU 方 案 对 2 m 温 度 和 2 m 相 对 湿 
度 的 模拟 效果 最 好 ,BouLac 方 案 对 10 m 风速 的 模拟 
效果 最 好 ,QNSE 方 案 对 PM 浓度 的 模拟 效果 最 好 ， 
整体 而 言 ,YSU 方案 对 民 勤 站 地 面 气象 要 素 的 模拟 
有 一 定 的 优势 ,QNSE 方案 对 这 些 要 素 的 模拟 效果 
相对 最 差 。 


参考 文献 (References ): 


[1] 张强. 大 气 边界 层 气象 学 研究 综述 四. 干旱 气象 , 2003, 21(3): 
74-78. [Zhang Qiang. Review of atmospheric boundary layer me- 
terology[J]. Arid Meteorology. 2003, 21(3): 74-78. ] 

[2] Banks R F, Tiana-Alsina J, Maria Baldasano J, et al. Sensitivity of 
boundary- layer variables to PBL schemes in the WRF model 
based on surface meteorological observations, lidar, and radio- 
sondes during the HygrA-CD campaign[J]. Atmospheric Research, 
2016, 176-177: 185-201. 

[3] Milovac J, Warrach- Sagi K, Behrendt A, et al. Investigation of 
PBL schemes combining the WRF model simulations with scan- 
ning water vapor differential absorption lidar measurements[J]. 
Journal of Geophysical Research: Atmospheres, 2016, 121(2) 
624-649. 

[4] Xie B, Fung J C H, Chan A, et al. Evaluation of nonlocal and local 
planetary boundary layer schemes in the WRF model[J]. Journal of 
Geophysical Research, 2012, 117(D12103). 

[5] IKEE, 刘 树 华 , 马 雁 军 , 等 . MYJ 和 YSU 方 案 对 WRF 边 界 层 气 
象 要 素 模 拟 的 影响 员 . 地 球 物理 学 报 , 2012, 55(7): 2239-2248. 
[Zhang Bihui, Liu Shuhua, Ma Yanjun, et al. The effect of MY] 


and YSU schemes on the simulation of boundary layer[J]. Chinese 


202103.00049v1 


chinaXiv 


T * 


Bl 


ChinaXiv 合 作 期 刊 


t 


E 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[5] 


[16] 


Journal of Geophysics, 2012, 55(7): 2239-2248. | 

黄 文彦 , 沈 新 勇 , 王 卫 国 , 等 . 边界 层 参 数 化 方案 对 边界 层 热力 
和 动力 结构 特征 影响 的 比较 町 . 地 球 物理 学 报 , 2014, 57(5): 
1399-1414. [Huang Wenyan, Shen Xinyong, Wang Weiguo, et al. 


Comparison of the thermal and dynamic structural characteristics 
in boundary layer with different boundary layer parameterizations 
[J]. Chinese Journal of Geophysics, 2014, 57(5): 1399-1414. ] 

EWE, EA, 赖 锡 柳 , 等 . WRF 模式 不 同 边界 层 参 数 化 方案 模 
拟 兰州 冬季 边界 层 高 度 的 研究 叫 . 高 原 气象 , 2017, 36(1): 162- 
172. | Wang Lixia Wang Ying, Lai Xiliu, et al. Study on the simu- 


lation of boundary layer height in Lanzhou in winter using WRF 
model with different boundary layer parameterization schemes|J]. 
Plateau Meteorology, 2017, 36(1): 162-172. ] 

HIR, OR R, 52448, 等 . 塔克拉玛干 沙漠 腹地 大 气 边界 层 参 
数 化 方案 的 模拟 评估 町 . 气象 科学 , 2018, 38(2): 157- 166. 


[Meng Lu, Zhao Tianliang, Yang Xinghua, et al. An assessment of 


atmospheric boundary layer schemes over the Taklimakan Desert 
hinterland[J]. Journal of the Meteorological Sciences, 2018, 38(2): 
157-166. | 

Gunwani P, Mohan M. Sensitivity of WRF model estimates to vari- 
ous PBL parameterizations in different climatic zones over India 
[J]. Atmospheric Research, 2017, 194(9): 43-65. 

Mohan M, Gupta M. Sensitivity of PBL parameterizations on PM i, 
and ozone simulation using chemical transport model WRF-Chem 
over a sub-tropical urban airshed in India[J]. Atmospheric Envi- 
ronment, 2018, 185: 53-63. 

Banks R F, Baldasano J M. Impact of WRF model PBL schemes 
on air quality simulations over Catalonia, Spain[J]. Science of The 
Total Environment, 2016, 572: 98-113. 

Sathyanadh A, Prabha T V, Balaji B, et al. Evaluation of WRF 
PBL parameterization schemes against direct observations during 
a dry event over the Ganges valley[J]. Atmospheric Research, 
2017, 193: 125-141. 

丁 成 慧 , 李 江南, 赵 杨 洁 , 等 . 边界 层 参 数 化 方案 对 南海 秋季 台 
JU TETUR" Q016) EET 09 57 LI]. 热带 气象 学 报 , 2018, 34(5): 
657-673. [Ding Chenghui, Li Jiangnan, Zhao Yangjie, et al. The 


influence of boundary layer parameterization schemes on autumn 
typhoon SARIKA(2016) in south China Sea[J]. Journal of Tropical 
Meteorology, 2018, 34(5): 657-673. ] 

温 晓 培 , 隆 雷 , 张 述 文 , 等 . 边界 层 参数 化 方案 对 台风 SANBA 9] 
生 阶 段 影 响 的 数值 模拟 研究 四. 热带 气象 学 报 , 2016, 32(3): 
346-357. [Wen Xiaopei , Long Xiao , Zhang Shuwen, et al. Param- 


eterization schemes on typhoon SANBA during its initial phase[J]. 
Journal of Tropical Meteorology, 2016, 32(3): 346-357. ] 

Liu J, Zhang F, Pu Z. Numerical simulation of the rapid intensifica- 
tion of hurricane Katrina(2005): Sensitivity to boundary layer pa- 
rameterization schemes[J]. Advances in Atmospheric Sciences, 
2017, 34(4): 482-496. 

崔 驰 潇 , 包 云 轩 , ORAS, 等 . 不 同 边界 层 参数 化 方案 对 江苏 地 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


p — WT ii 29. iE Fe BRE TELE WS] [3]. 大 气 科学 , 2018, 42(6): 
1344-1362. [Cui Chixiao, Bao Yunxuan, Yuan Chengsong, et al. 
Influence of different boundary layer parameterization schemes on 
the simulation of an advection fog process in Jiangsu[J]. Chinese 
Journal of Atmospheric Sciences, 2018, 42 (6): 1344-1362. ] 
Mylonas M P, Nastos P T, Matsangouras I T. PBL parameterization 
schemes sensitivity analysis on WRF modeling of a tornadic event 
environment in Skala Lakonia in September 2015[J]. Atmospheric 
Research, 2018, 208: 116-131. 

18 2508, 李 牙 清 , EX, 等 . WRE 模式 多 种 边界 层 参数 化 方案 
对 四 川 盆 地 不 同 量 级 降水 影响 的 数值 试验 中. 大 气 科学 ， 
2016, 40(2): 371-389. [Gao Duming, Li Yueqing, Jiang Xingwen, 
et al. Influence of planetary boundary layer parameterization 
schemes on the prediction of rainfall with different magnitudes in 
the Sichuan Basin using the WRF model[J]. Chinese Journal of At- 
mospheric Sciences, 2016, 40(2): 371-389. ] 

ZE3E, AS, 周立波 , 等 . WRF 模式 中 边界 层 参 数 化 方案 在 藏 东 
南 复杂 下 垫 面 适用 性 研究 川 . 高 原 气 象 , 2017, 36(2): 340-357. 


[Li Fei, Zou Han, Zhou Libo, et al. Study of boundary layer param- 


eterization schemes ' applicability of WRF model over complex un- 
derlying surfaces in Southeast Tibet[J|. Plateau Meteorology, 
2017, 36(2): 340-357. ] 

Wang J K, Zhang B H, Zhang H D, et al. The impacts of three PBL 
schemes on the simulation of boundary layer meteorological fac- 
tors and air pollutants in Northern China[JJ. IOP Conference Se- 
ries: Earth and Environmental Science, 2019, 227(5): 98-113. 
ACER, 王 金 艳 , 气 继 勇 , 等 . 起 沙 方案 对 西北 地 区 沙 尘 过 程 模 
拟 的 影响 四 . 环境 保护 科学 , 2018, 44(4): 69-76. [Liu Xiaoran, 


Wang Jinyan, Qiu Jiyong, et al. Impact of dust emission schemes 


on the simulation of dust storms in the Northwest region[J]. Envi- 
ronmental Protection Science, 2018, 44(4): 69-76. | 

TEAN. 基于 超大 尺度 结构 模 化 的 沙尘暴 及 沙 墙 形态 数值 研究 
[D]. 兰州 : 兰州 大 学 , 2018. [Huang Haojie. Numerical Studies on 


Sand-Dust Storms and Sand Wall Morphology Based on the Model- 
ing of the Very-Large-Scale Motions[D]. Lanzhou: Lanzhou Univer- 
sity, 2018. ] 

赵建华 , 张 峰 , RE SE. GIU CHE T A FATE EUER]. 
干旱 区 研究 , 2019, 36(6): 1419-1430. [Zhao Jianhua, Zhang 


Feng, Liang Yun, et al. Research progress on turbulent coherent 
structure in atmospheric boundary layer[J]. Arid Zone Research, 
2019, 36(6): 1419-1430. ] 

Hong S Y. A new vertical diffusion package with an explicit treat- 
ment of entrainment processes[J]. Monthly Weather Review, 2006, 
134(9): 2318-2341. 

Janjié Z I. The Step-Mountain Eta Coordinate Model: Further devel- 
opments of the convection, viscous sublayer, and turbulence closure 
schemes[J]. Monthly Weather Review, 1994, 122(5): 927—945. 
Nakanishi M, Niino H. An improved Mellor-Yamada Level-3 Mod- 


el: Its numerical stability and application to a regional prediction 


202103.00049v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


TRO 全 等 :边界 层 参 数 化 方案 对 一 次 西北 地 区 沙 侍 天 气 过 程 影响 的 数值 模拟 研究 


of advection fog[J]. Boundary-Layer Meteorology, 2006, 119(2): [J]. 兰州 大 学 学 报 ( 自 然 科 学 版 ), 2012, 48(1): 69-74. [Long 
397-407. Xiao, Zhao Jianhua, Huang Jianping. Observation and analysis of 
[27] Sukoriansky S, Galperin B, Perov V. Application of a new spectral *07. 3° severe sandstorm[J]. Journal of Lanzhou University (Natu- 
theory of stably stratified turbulence to the atmospheric boundary ral Sciences), 2012, 48(1): 69-74. ] 
layer over sea ice[J]. Boundary- Layer Meteorology, 2005, 117(2): [31] HA, EE TS, FU, 等 . 民 勤 一 次 沙尘暴 天 气 过 程 的 近 地 层 气 象 
231-251. HRE KERE. 干旱 区 人 研究, 2018, 35(6): 1352-1362. [Wei 
[28] Bougeault P, Lacarrre P. Parametrization of orography-induced tur- Qian, Long Xiao, Tian Chang, et al. Multiscale meteorological 
bulence in a meso-beta model[J]. Monthly Weather Review, 1989, characteristics during a sandstorm in Minqin[JJ. Arid Zone Re- 
117: 1872-1890. search, 2018, 35(6): 1352-1362. n 
[29] Nabavi S O, Haimberger L, Samimi C. Sensitivity of WRF- chem [32] Sfi. 一 次 沙 尘 天 气 过 程 发 展演 变 的 多 尺度 特征 研究 [D]. 兰 
predictions to dust source function specification in West Asia[J]. 州 : 兰州 大 学 , 2020. [Wei Qian. Studies on Multiscale Character- 
Aeolian Research, 2017, 24: 115-131. istics during the Evolution of A Dust Event[D]. Lanzhou: Lanzhou 
[30] HET, 赵建华 , 黄建平 .“07. 3? 强 沙尘暴 天 气 过 程 的 观测 分 析 University, 2020. ] 


Impact of boundary layer parameterization schemes on the 
simulation of a dust event over Northwest China 


WEI Qian, LONG Xiao, ZHAO Jianhua, HAN Zifei, WANG Siyi' 
(1. College of Atmospheric Sciences, Lanzhou University, Lanzhou 730000, Gansu, China; 
2. Institute of Arid Meteorology, China Meteorological Administration, Lanzhou 730020, Gansu, China) 


Abstract: In this study, WRF-Chem(version 3.4) model was used to compare the performance of different 
Planetary Boundary Layer (PBL) parameterization namely, the YSU(Yonsei University), MYJ(Mellor-Yamada- 
Janjic), QNSE(Quasi- Normal Scale Elimination, MYNN2.5(Mellor- Yamada- Nakanisfii and Niino 2.5) and 
BouLac PBL schemes, over the dust event in northwest china on 27 March 2007. Surface observations were used 
for comparisons and evaluating model performance for meteorological variables. It is shown that simulations 
with the five PBL schemes can successfully reproduce the evolution of the dust event. The YSU and BouLac 
schemes produced higher surface friction velocity, 10 m wind speed, 2 m air temperature and surface PM 
concentration and lower 2-mrelative humidity, thus simulating stronger weather processes than those of the MY]J, 
QNSE and MYNN2.5 schemes. These results indicate that different boundary layer schemes affect the dust 
emission flux and PMi concentration through different simulation effects of friction velocity. The dust emission 
flux and PM, concentration tend to increase with higher friction velocity. Therefore, the dust event was enhanced 
due to the high friction velocity and the characteristics of high temperature as well as lower humidity in near- 
surface layer in the afternoon. As a result, the simulated dust event with the BouLac scheme was the strongest 
while the weakest by the QNSE scheme. Observations from Minqin meteorological station are used to validate 
the simulated results over Mingqin region. Statistical analysis of the five simulations shows that the QNSE scheme 
simulated better PMio concentration, the BouLac scheme performs well for 10 m wind speed and the YSU scheme 
resulted in the best model performance for simulating air temperature and relative humidity at 2 m. Overall, the 
YSU scheme was concluded as the best PBL scheme for the dust storm and the QNSE scheme was the worst one. 


Keywords: WRF-Chem; sand and dust; planetary boundary layer; numerical simulation 


